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Abstract 

We present a systematic derivation of the Biot-Savart equation from the Nonlinear Schrodinger 
equation, in the limit when the curvature radius of vortex lines and the inter-vortex distance are 
much greater than the vortex healing length, or core radius. We derive the Biot-Savart equations 
in Hamiltonian form with Hamiltonian expressed in terms of vortex lines, 

with cut-off length C ~ 0.3416293/^^, where po is the background condensate density far from 
the vortex lines and k is the quantum of circulation. 
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I. INTRODUCTION 


Nonlinear Schrodinger (NLS) equation Q and Biot-Savart equation Q are the two most 
popular models for describing superfluid dynamics and turbulence at very low temperature 
m- This includes flows in superfluid ^He and ^He, as well as atomic Bose-Einstein con¬ 
densates of alkali gases. It is generally believed that the Biot-Savart equation describes a 
subset of slow (subsonic) motions of the more general NLS equation. This view is based on 
the formal ideal-fluid analogy in the NLS model arising from the Madelung transformation, 
as will be explained below in Section [ITB However, to date there has been no rigorous jus- 
tihcation for such a correspondence between the Biot-Savart and the NLS equations. The 
problem is that the zero vortex line radius limit is ill dehned for the ideal fluids leading to a 
divergence of the Biot-Savart integral. To tackle this problem, it is customary to introduce 
a phenomenological cut-off regularisation of the Biot-Savart integral at a scale that corre¬ 
sponds roughly to the vortex line radius. All efforts to justify this approach rigorously for 
the ideal fluid dynamics have failed due to a great freedom in possible realisations of the 
vortex prohles and their temporal variability due to vortex stretching. 

Fortunately, the fluid dynamics equations arising from the NLS equation contain an extra 


term - the so called quantum pressure (see Section IIB). It is the quantum pressure term 
that makes the quantum vortex core “rigid”, i.e. having a hxed universal prohle. This fact 
makes it possible to rigorously derive the Biot-Savart equation from the NLS equation. The 
present paper is devoted to such a derivation. 


II. NONLINEAR SCHRODINGER EQUATION 

Consider a model described by the 3D defocusing Nonlinear Schrodinger (NLS) equation 

PHD], 

ihdr'^ + + FT - = 0 , 

where T is a complex scalar function of a 3D physical space coordinate X and time T, h is the 
reduced Planck’s constant, M and E are the mass and single-particle energy of the bosons, 
and Vq is the strength of the interaction potential between them. Dehne dimensionless 
variables as follows: x = t = and ip = The resulting equation is: 

+ = 0 . ( 1 ) 
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A. Hamiltonian formulation 


The NTS equation Q can be written in Hamiltonian form, 

6H 


idtilJ = 


S'ljj* ’ 


with Hamiltonian 


H = 


|VV’(x,t)|2 + - (|V’(x,t)P - 1)' 


(ix, 


( 2 ) 


( 3 ) 


which represents the conserved energy of the system. 


B. Madelung transformation and fluid framework 


The Madelung transformation P, [TH |T2] maps the complex scalar field '0(x, t) to two 
real scalar helds p(x, f) and 0(x, t) following the relation ■0 = By plugging this 

substitution into the defocusing NLS equation ([^ and separating the real and imaginary 
parts, one obtains the set of equations 

|^ + V.(p 2 V^) = 0 , 

^ + (v.#.r + p-i-^ = o. 

It is then straightforward to observe that by setting the vector held u(x, t) = 2V(j){'x, t) the 
hrst equation results in a continuity equation for the density held p and the second equation 
a conservation of the momentum associated with the velocity held u: 


| + V.(pu)=0, 
| + (u.V,u = 




p \ \fp 

Thus, we have obtained equations of an inviscid polytropic gas with adiabatic index 7 = 2 
(pressure p = p^), but with an extra term, called quantum pressure, appearing at the last 
term in the second equation. The quantum pressure term is negligible if the characteristic 
scale of motion is much greater than the healing length 

Let us re-express Hamiltonian ([^ in terms of the huid variables: 


H^\ 


(p - 1)^ + 2|Vv^P 


dx. 


( 4 ) 


With the exception of the last term, we see a formal coincidence with the standard expres¬ 
sions for the total energy (up to factor 1 / 2 ) of a compressible huid. 
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III. QUANTISED VORTICES 


Even if the velocity field u is irrotational, vortices may appear in the system. These 
are lines of singular vorticity around which the real phase held changes by a multiple of 2tt 
assuring that the complex wave held '0 stays always single-valued. One can measure the 
circulation along a closed curve C around one of these vortex lines 

r = ® u - dl = 2 (p V0 ■ dl = 2 Acj) = dvrn, with n G Z . 

Jc Jc 

The circulation is zero for all contours embracing the regions where the phase is a well- 
dehned diherentiable function. However, at points where "0 = 0, the phase is undehned and 
the circulation along contours embracing such points is not zero. The circulation may take 
only discrete values and the singularities in the vorticity held are called quantum vortices. 
Moreover, vortices with n > 2 are unstable: a general smooth change in held 'ijj leads to 
splitting of such “multi-charge” vortices into a set of elementary vortices with n = 1. Because 
of the quantised circulation, together with the fact that the huid has no viscosity, the NTS 
model has been used to qualitatively describe superhuids. 

Let us hrst consider a single straight vortex solution found by Pitaevskii [10]. For sim¬ 
plicity, let us consider a situation where the density far from vortices asymptotes to po = 1- 
The more general case when the density tends to a diherent constant po > 0 can then be 
obtained by a simple rescaling '0(x, t) —)• y/po4’i\fPo^jt)- Let us impose that a phase shift 
of +271 exists around the origin for the phase held 0 and look for a stationary solution. 
Going into polar coordinates 

{ X = r cos 9 
y = r sin 6 

we consider a solution with 0(r, 9) = 9. We will also impose an axial symmetry on the 
density held, so that the solution has the form 


iiv{r,9) = R{r) e 


ie 


( 5 ) 


Plugging this into the NLS equation 0 we get the ordinary diherential equation 

d^R 

-I- 

r dr 


1 dR 1 „ . ^ 

— — - — = 

T* /Y T* rr^ ^ 


( 6 ) 


df2 j. 

supplemented with the boundary conditions i?(0) = 0, R{oo) = 1. One may obtain ana¬ 
lytically the behaviour of the function R{r) in the limits of very small and very large radius. 
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or look for a numerical solution using a shooting method. Another method is to use a Fade 
approximation [13]. For the present paper, it turns out, the computation of the cut-off 
length (see Section 0 requires knowledge of a very accurate solution for R{r) over a wide 
range of values of r. This rules out Fade approximations (they have signihcant errors of up 
to 3.5%) so we are forced to work with high-accuracy numerical solutions and asymptotic 
solutions at large r. Figure shows a log-linear plot of R{r) obtained by our method as 
explained in Appendix]^ The vortex core has size of the order of the healing length tending 
rapidly to a constant value far from the origin. The held V’ is smooth everywhere and tends 
to zero at the origin which can be thought of as a localised topological phase defect. 

R(r) 



FIG. 1. (Color online) Log-linear plot of the vortex profile R{r), obtained as a 
nation between a highly-accurate numerical solution and an asymptotic solution, 
and Supplemental Material m for details. 


piecewise combi- 
See Appendix |A| 


IV. VORTEX TANGLE AND BIOT-SAVART MODEL 

General vortex structure in strong 3D NLS turbulence may be very complex and irregular: 
it is usually referred as vortex tangle. In a wider context, vortex tangle represents a typical 
realisation of superhuid turbulence at zero temperature, i.e. in liquid Helium 4 or 3. 

When the distance between the vortex lines and their curvature radii are much greater 
than the healing length the hrst term in the integral of energy (|^, the so-called ki¬ 
netic energy, is dominant over the second and the third contributions, the internal and 
the quantum-pressure energies respectively. This is because (provided there is no ambient 
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sound) the density field deviations from the background density po (in our case po = 1) are 
bound within small vortex cores of radius ~ whereas the velocity field u produced by the 
vortices is de-localised. Thus, the leading order in the Hamiltonian is 


H f pu^ dx, 


(7) 


which is a Hamiltonian of an incompressible fluid (up to factor 1/2), but with density 
depletions in the vortex cores which are rigidly fixed to have transverse density distributions 


of the isolated 2D vortex considered in Section HI (Note for our purposes we will also need 
to find the main sub-leading order in id). 

This suggests that the vortex tangle can be modelled by the Biot-Savart equation, which 
has been used the context of the vortex dynamics in incompressible fluids starting with the 
work of Da Rios 1151: 


_ ds K 

* dt dvr 


(s' — s) X ds' 

Is — s'13 


( 8 ) 


where s = s{(,t), where G M is a Lagrangian label parametrising the positions of the 
vortex line elements. Here, k is the quantum of circulation; in the non-dimensional NLS 
model considered in this paper we have k = dvr. The integral is taken along all of the 
vortex lines, including the one to which the considered vortex element belongs. To prevent 
logarithmic singularity at s' —s, the integration has a cut-off at the “vortex radius” scale, 
< |s — s'l, whose value, based on the common sense physical grounds, must be close to 
the healing length Such a cut-off has been previously introduced phenomenologically, but 
later in the present paper we will provide a rigorous justification for this model and will find 
the value of numerically. 

Equation (|^ can be obtained from the least-action principle 

6S= f 6Cdt = 0 


with Lagrangian 


and Hamiltonian 


D = ^ J Sf (s X ds) —T-L 


n = 


K 


^ f ds' ■ ds'‘ 


Ibvr / Is' — s" 


( 9 ) 


Variation of action gives: 


6S = 


dt (5s ■ 


K 5H 


dC = o. 
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Considering the fact that the parametrisation ( is arbitrary, the expression in the square 
brackets must be zero for any vector tangential to s. The respective Hamiltonian equation 
is: 

-s,xs, = --. (10) 

Taking into account that 


5H ^2 

f (sc- (s-s'))ds'- ( 

s^ ■ ds'){s — s') 

/■ Sc X [(s' - s 

) X ds'] 

Ss Svr j 


s — s' 

|3 StT j 


s — s' 

3 


and undoing the cross product, we get the Biot-Savart equation (|^. 


V. DERIVATION OF THE BIOT-SAVART MODEL FROM THE NLS EQUATION 


In spite of great popularity of the Biot-Savart model, it has not yet been rigorously 
obtained or justihed in the form formulated here, i.e. with an integration cut-off. In the 
ideal fluids context, the difficulty is related to the vortex stretching effect, which causes 
changes in the vortex core radius. The situation may be easier in the NLS model, since the 
vortex core shape is hxed. However, previous attempts to derive the vortex line evolution 
in the NLS model have led to signihcantly more complex equations than equation ([^ (see 
e.g. [I9]). 


The present paper is devoted to such a derivation of the Biot-Savart model within the 
NLS model. For this one has to (A) rewrite the NLS equation in the vortex hlament form 


(10), and (B) express the Hamiltonian H in terms of the vortex line configuration 
Part (A) was done in [8] in more general settings (forced Ginzburg-Landau equation), and 
we reproduce it for our NLS case in the Appendix On the other hand, part (B) has 
not been done before to the extent that the cut-off in the Biot-Savart equation ([^ and the 
Hamiltonian ([^ would be rigorously justihed and calculated, and this will be the main goal 
of the derivations that follow. 

The process of hnding the Hamiltonian in terms of the vortex line is done in six steps. 
The result will be 


H = 


K 

Ibvr 


|s-s'|>5. 


(is ■ ds' 
Is — s'l 


( 11 ) 


where 


0.3416293 ± 10 The value of the cut-off will be shown to satisfy the following 
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analytical formula: 




r dr — lim 

r—^oo 



[Rjrr 



where R{r) is the vortex prohle, solution of equation ([^. 


( 12 ) 


Step 1: Obtaining the vortex-line Hamiltonian. This was already done by writing 
(|^. We will assume no sound is present or excited by the moving vortices, which itself 
requires a justihcation if one wants to be fully rigorous. However, we expect this to be a 
difficult task akin to proving the validity of the balanced (wave-free) geophysical motions. 
Physically, we could say that sound is not generated when the motions of the vortices are 
strongly subsonic, which is true when they are separated by the distances much greater than 
their core radii. 

Thus, here we will simply postulate Hamiltonian (|^ and assume that the density field 
rigidly follows the vortex cores, so that locally near the cores we have the Pitaevskii vortex 
prohle. In this way, the last two terms in this Hamiltonian will be sub-dominant (after 
subtracting the constant contribution), but still important because they affect the value of 
the cut-off length 


Step 2: Reformulating the problem in the form of a flow with constant density 
in the leading order. We do this by introducing a new “velocity” held 

(13) 

SO that the Hamiltonian is 

H = Hk + Ho, HK = \jv^d^, Ho = ^J[{p-lY + 2\V^\^]d^, (14) 

where Hk formally looks like the kinetic energy of incompressible fluid with density equal 
to one. The terms in Hq depend on the vortex profile and contribute directly to the value 
of the cut-off length 

For an infinite straight (Pitaevskii) vortex, we have: 










where 6 is the unit vector in the azimuthal direction and rj_ is the distance to the vortex 
line. 

Taking into account that that a/ p/2 —)■ ar±, a = const, when the distance to the vortex 
rj_ —>■ 0, and that in the vortex u = 2/r±, we have 

V —)■ 2a, as r±_ —)■ 0. (15) 


Same relation holds for a more general vortex line when its curvature radius is much greater 
than the healing length. 

We compute Hq explicitly in the case of the inhnite straight vortex. This will be used 
later on. In terms of the vortex prohle R{r) we obtain, per unit length of vortex line. 


_ k" 
to'*"' 


where k is the quantum of circulation, k = Att, and 


/^o = 


'0 




r dr. 


These integrals can be split as follows. First 

poo 

/ (\R{ryf‘— lY r dr = 1. 

Jo 

Second, we obtain numerically 


dR{r] 

dr 


r dr ^ 0.279090913. 


Therefore 


(16) 


(17) 


lio ~ 0.779090913 


Step 3: Rewriting the Hamiltonian (14) in such a way that the integrand func¬ 


tions are localised within the vortex cores. This can be done at the expense of 
increasing dimensionality of the integral via performing an integration by parts followed by 
applying the Helmholtz theorem: 


Hk = — 
Svr 


1 /■ u;(x, t) ■ u;(x', t) + 7 (x, t) 7 (x', t) 


X — X' 


d:sidx', 


(18) 
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where <-j(x, t) is the “u-vorticity” field: <-j(x, t) = V x v(x, t), and 7 (x, t) is the n-divergence 
field: 7 (x, t) = V ■ v(x, t). 

To obtain this formula, we have had to use Helmholtz theorem on a domain that does 
not include the vortex lines, so that v(x, t) satisfies the regularity conditions of the theorem. 
Now, it is possible to make this domain as close to the vortex lines as we like and show 
that the contribution of the remaining small vicinity of the vortex is vanishing because the 


singularity of v(x, f) is mild enough. So the integration domain in equation (18) is the full 


From here on, we will discard the contribution from the divergence terms 7 (x, t), because 
they are small if the vortex line’s curvature radius is much greater than the healing length. 


From (15) we have asymptotics for the vorticity near the vortex centre: 


1 / 2a 

u = - Or, [r±v) —)• —, as r± —)■ 0. 
r± r± 


(19) 


The vorticity field is strongly localised within the healing length from the vortex centre and 
rapidly decays at large distances from the centre. 

We see that such a n-vorticity is no longer a delta-function distributed on the vortex 
centre, as it was the case for the original w-vorticity. It is still singular at r± = 0, but 
the singularity is mild, and it leads to a non-divergent H —without a cut-off. The latter 
statement is true because the integral in H is identical to the original convergent integral ([^. 


Step 4: Splitting the Hamiltonian in terms of local and far regions. This is the 
key step. We introduce an intermediate length scale a which is in between of the healing 
length ^ and the curvature of the vortex line t. ^ a <^ £. We split the integrals in Hk as 




where 


and 


H> = — 

Svr 


H^ = — 


Stt 


|x—x'|>a 


' |x—x'|<a 





^(x. 

X — x' 



X — x' 



SO i/> and iJ< are the integrals over the domains |x — x'| > a and |x — x'| < a respectively. 
Notice that the term Hq contributes to H with a term that is more in the spirit of if< than 
Hy because the integrand in Hq is supported in the vicinity of the vortex lines. 
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Writing i/> in the Biot-Savart form Q comes cheaply because one can replace |x — x'| 
in the denominator with |s — s'|, which would be valid up to 0{^/a) corrections. Then the 
transverse to the vortex line directions can be integrated out independently for the primed 
and the un-primed variables, and in the end we have: 

^ ds ■ ds' 


= 


K 

Ibvr J^s—s'\>a I® ® I 


( 20 ) 


where k is the quantum of circulation, k = Ati. The pre-factor here comes from the fact that 
the transverse integration of tc;(x) gives the circulation of v(x), i.e. 2^/^7r. 

To consider the if< integral, one can think of the vortex line as locally straight and 
use local coordinates. Let us consider two planes A and A' transverse to the vortex line 
and passing through the points x and x' respectively. Let us introduce the change of the 
integration variables (x, x') —(sy, r_i_, sy', r'_j_), where sy and sy' are the lengths along the 
vortex line to its intersections with planes A and W; and r'j^ are the local 2D Cartesian 
coordinates within A and A' respectively. Taking into account that the contributions to the 
integral are limited to local regions, |x — x'| < a, the change of variables amounts to a shift 
and rotation, i.e. its Jacobian is one, up to 0(a/£) corrections. The integration region is 
then |sy — Sm| < a, and we have: 






T||-sj|l<a y^(sj| - Sy)‘ 


dswdr I dsudr'i 


+ (r', 




After integrating out r_i_ and r'j_, the remaining integrand will be a function of s = |sy — sy| 
only (this s should not be confused with |s|). Also, up to 0{a/1) corrections integral Hq can 


be replaced with its value obtained for the straight line, equation (16). Combining Hq and 
J/<, we get 


TJ ZJ ^ 

Lz< -l- Hq — -^^Hq -\- 


Stt 


IGtt 


b||-dll<a 


/(|sy — Sy I) ds\\ds'u = 


k^L 

Stt 


ho + / /(s) ds , (21) 


where /xq is a constant dehned in equation (17), L is the total length of the vortex hlament, 
and 

u){r±)u{r'j 


/(«) = 


idr I dr'i 


To obtain the latter expression we performed a linear transformation from (sy,Sy) to 
(sy — Sy, Sy + Sy), with Jacobian equal to 1/2. 
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Step 5: Defining the effective cut-off analytically and finding the effective Hamil¬ 
tonian. At this point it is useful to define the function 


F{a) = /io + / /(s) ds. 

Jo 

At large distances, 1 ~ ^ <C s ^ a, we have /(s) ~ 1/s so that 

F{a) =fio+ /(s) ds ^ Ina + fio + C = ln(a/^*). 
Jo 


( 22 ) 


(23) 


This can be considered a definition of the effective cut-off length in the Biot-Savart 
formulation, which is expected to be different from the standard expression .^ = 1/= 1 
by an order-of-one factor. The task of finding numerically is left to the next step. Notice 


that at s —)• 0, function /(s) = o(l/s), so the integral in (23) converges and T"(0) = /iq. In 
fact, /(O) is finite (~ 1.122 by numerical estimation). 

To obtain a closed formula for the full Hamiltonian, define the function 


g{s) = s/(s) = ^ 




, _idridr'i , 

a/s 2 (r'_L - r_L)2 

which represents a “smooth cut-off”; (/(O) = 0, g{oo) = 1 (see figure [^. 


(24) 


From equations (21), (22) we rewrite 

k^L 




H. = 


Stt 


F(a) = 


k‘L 

Stt 


ds 


K 


ds ■ ds' 




Ja>\s—s'\>^t I® ® I 


Adding this expression to Hy given in Eq.(20) we have: 

5 f(|s — s'l) ds ■ ds' 


H = —Aio 


K 


ds ■ ds' 


s — s' 


16vr |s s I 


Stt IGtt J 

which is Eq.(pT| we were aiming to derive. This is an asymptotically exact result valid when 
the vortex line’s curvature radius is much greater than the healing length. 


Step 6: Finding the cut-off length numerically. To compute 5 '(s), we rewrite (24) 
in polar coordinates: 


d(«) = f 


r ±r' ±uj{r [)u{r'^) 


dOdr I dr'i 


47r J y/s2 _|_ _|_ ^2 _ 2r_i_r(|_ cos 9 

The angular integration can be performed analytically, leading to an expression for g{s) that 
is easier to handle numerically. Dropping the ‘perp’ symbols we get: 


(*OQ /*00 


K 


9{s) = - / 

TT Jo Jo 


Arr' 


r oj{r) r'oj{r') — ^ ^ drdr', 

+ r'Y + 
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where K{m) is the complete elliptic integral of the first kind. Eqs. (13) and (19) lead to 
relation ru{r) = a/2 , where R{r) is the solution of the ODE ([^ (in terms of the vortex 

density we have R{r) = p{r)). Thus, we get: 


=V r 

^ io 


dR{r) dR{r-') ^ l^(r+r') 2 +s 2 


Arr 


dr 


dr' ^(r+ r')2 + s2 


drdr'. 


(25) 


We will compute this integral numerically using a very accurate numerical solution of 
equation ([^ for R{r), combined with an asymptotic solution valid for large r, so that the 
error incurred in solving equation (§ is kept uniformly below 5 x 10 (see Appendix 
A). Replacing the corresponding expression for into (25), the resulting 2D integral is 


computed numerically for several values of s. The partial result is plotted in hgure 

To obtain the effective cut-off length we need to integrate numerically /(s) (= g{s)/s) 
over the range s G (0,a), thus obtaining F{a) (= po + Jq f{^)ds), and compare the result 


with Eq.(23). The numerical computation of po is done accurately using the prohle R{r) 
and the result is po ~ 0.779090913. Figure shows the results of numerically integrating 
/(s) in segments, using Mathematica's global adaptive method with accuracy and precision 
goals of 10“^° each. The blue dots represent the values of f{s)ds obtained numerically. 
The red straight oblique line corresponds to the £t // f{s)ds ~ log(a) -|- C which is good i 


m 


the asymptotic regime a S> 1 (we went up to a > 10® to verify the asymptote). By looking 
at figure the value exp(—C) is by definition the value of a at which the horizontal line 
and the red straight oblique line intersect. Combined with the numerical value of po, the fit 
gives the following accurate estimation for the cut-off length: 


= exp(-po -C)= 0.3416293 ± 10 


-7 


( 26 ) 


This result can be immediately generalised to the case po / 1: 


^ 0.3416293/7)5^. 


(27) 


Finally, in Appendixj^we derive the analytical formula (12) that allows us to bypass the 


numerical fitting procedure. Using the accurate numerical solution for R{r) on formula (12) 


yields the same value for 7 as in equation (27) above 
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g(s) 

1 . 0 - 



FIG. 3. (Color online) Blue filled circles: Numerical computation (in segments) of the integral 
Jq f{s)ds. The red oblique straight line corresponds to the ht Jq fis)ds ^ ln(a) + C. The green 
vertical dot-dashed line corresponds to a = exp(—C). 


VI. CONCLUSIONS 


In this paper, for the hrst time, we have rigorously derived the Biot-Savart equation ( pT| 
with cut-off (12) from the 3D defocusing NLS equation ([^. For this we assumed that the 
vortex line curvature and the inter-vortex distance are much greater than the healing length. 
This setup corresponds to a very subsonic motion of the vortex lines, such that the density 
far from the vortices tends to a constant value po- Correspondingly, we assumed that the 
core is rigid and has a hxed angle-independent prohle. This is justihed from the fact that all 
angle-dependent modes quickly dissipate via sound radiation |20]. Non-circular corrections 
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to the core shape arising from curvature (e.g. for a vortex ring) are small and can be ignored 
in our derivations. 

We have found an accurate numerical value for the cut-off length, 0.3416293/y^. 
This value agrees with the common sense suggestion that it must be of the order of the 
healing length ^ = 1/y/po. However, the numerical coefficient shows an approximately 
threefold difference between these two quantities, which may introduce signihcant corrections 
to previous estimates of, e.g., the critical velocity for the vortex formation in a superfluid 
flow past an obstacle, critical distance for the vortex line reconnection process, speed of 
vortex rings, etc. 

It is important to understand that the NTS model provides us with a signihcant advantage 
with respect to the ideal huid model described by the Euler equation: the vortex prohle is 
hxed in the former but not in the latter. As a result, the Biot-Savart equation would be 
much harder to derive from the Euler equation, and it would be unrealistic to expect that 
the result would be in a simple form as in Eq.(pT| with a hxed cut-oh length. However, 
it would be interesting to explore possible justihcations of the Biot-Savart model extension 
with a variable cut-oh length physically corresponding to variability of the vortex radius due 
to the vortex stretching process. 

In the future, it would be interesting to conduct a comparative numerical study of the 
Biot-Savart and the NLS models for various physical examples, including propagating and 
colliding vortex rings, Kelvin waves on vortex lines and reconnection of vortex lines. In 
particular, it would be interesting to see at what stage and under which conditions the 
Biot-Savart approach fails due to, e.g., close approach of diherent vortex lines to each other, 
appearance of regions with high curvature or loss of energy to acoustic waves. 
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Appendix A: Accurate solution of vortex profile R{r) 


We implemented a very accurate shooting method to obtain R{r) numerically between 
r = tq = exp(—50) ~ 10“^^ and r = 20. The method uses Mathematica's stiffness switching 
solver with working precision of 22 digits and accuracy and precision goals of 10 “^^ each. 
As a result, the absolute point-wise errors in solving equation ([^ is kept below 10 “^^ for all 
r in the range 10“^^ < r < 10“^, and below 10“^^ in the range 10“^ < r < 20. The optimal 
value of derivative at ro is found to be 

= 0.5831894958603292791791737551. 

r=ro 

It turns out that this numerical solution cannot be continued much beyond r = 20 without 
significant loss of accuracy. In order to continue the solution to higher values of r, we devised 
an asymptotic method. This method consists of the transformation 

R{r) = exp(—Z(l/r^)), 

which leads to the ODE 

1 _ p- 2 .Z(g) 

- Aq’^Z'iq? + 4qZ'(q) -+ 1 = 0. 

q 

This ODE is solved near g = 0 in power series 

N 

Z{q) = cjq^ , 

j=i 

leading to an asymptotic solution that does not converge, but can be truncated so that the 
error in solving equation ([^ is kept as small as desired for r large enough. In practice we 
found that A = 13 is a good compromise. The coefficients ci, C 2 ,..., C 13 are, explicitly: 

1 5 32 1589 64981 989939 168211250 38006710085 5510235057787 
2’ 4’ y ’ ~ 8 ~’ 10 ’ 3 ’ 7 ’ 16 ’ 18 ’ 

199454257136329 110192683498843556 14600012068277445755 9139380150115822460510 
4 ’ n ’ 6 ’ 13 ■ 

Next, we patch this asymptotic series with the previous numerical solution by finding the 

intersection of the two functions. They intersect at 


dR 

dr 


r-trans = 15.5752238612175596563309, 
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which defines a piecewise continuous solution at the transition point utrans- The jump in the 
derivative of this piecewise function is found to be reasonably small. In absolute terms, 


dR 


dR 


dr 

^ ^trans 

dr 

^ ^trans 


< 2 X 10■^^ 


while in relative terms, 


dR 


dR 


dr 

^ ^trails 

dr 

^ ^trans 


< 5 X 10-^2 


dR 


dr 


r=r 


+ 

trans 


These errors are mostly due to the asymptotic solution’s error in solving equation ([^, which 
is of order 2 x at the transition point r = rtrans and rapidly decreases below 10“^^ for 

r > 21. 

Details of the implementation of this accurate piecewise solution are found in the Sup¬ 
plemental Material na. 


Appendix B: Hamiltonian eqnations in terms of vortex lines 


Below, we reproduce derivations of Nemirovskii [8] simplifying them as appropriate for 
the NLS model Q. Consistently with what we assumed when deriving H, we will consider 
highly subsonic motions of the vortex lines which occur when their curvature radius and 
mutual separations remain much greater than the healing length In this case one can 
neglect the acoustic waves and assume that function ipiyiR) is fully determined by the 
vortex line configuration 

= V’(x|s(C,t)). 


In particular, the time dependence in appears only implicitly via s{(,t) so that 




dC'. 


(Bl) 


5s(C',t) 

Let us multiply equation (|^ by /5s{(Q,t) add its complex conjugate and 

integrate over the 3D physical space: 


i / dt'ijj{x,t) 


dip* 


(5s(Co,t) 


dx -|- c.c. = 


5H dip* dH dip 


+ 


dx = 


dH 


(B2) 


is(Ca,t) ' SiiSs{(a,t)J Ss{(a,t)' 
Dominant contribution to the integral on the LHS of this equation comes from a small 
vicinity of the vortex line, where locally ip behaves as the Pitaevskii vortex ([^, i.e. 

^(x|s(C, t)) ^ 1 pv{x±) = '0^(s(C±, t) - x). 
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where corresponds to the point on the vortex line which is closest to x. Then we can 
write 


(5V>(x|s(C,t)) 

6s{C,t) 


= Vx^.(x'Jh(C'-Ca¬ 


using this in the LHS of (B2) and its conjugate version in (Bl), and integrating over 


dx\\ = |s^(Co, t)|(iC± and over (' using the delta-functions, the LHS of (B2) becomes: 


i|s^(Co,t)| / (st(Co,t) • V±V'^(x_L))V±^*(x_L)dx_L-hc.c. 


Using the vector triple product formula, we have: 


i (St(Co,t) ■ V±V't;(x±))V±^/’*(x_L) C.c. =iSt{Co,t) X (V±'0i^(x_L) X V_L'0i^(x_L)) 


Substituting ([^, we have 


V±V^*(x_l) X V±V’t;(x±) = t. 


x± 

sc(Co.b 


where the prime here means derivative, x±_ = |xj_| and t = is a unit vector tangential 

to the vortex line. We have taken into account that V_l6 * = 9/x±^ where 0 is a unit vector 
in the azimuthal direction. 


Thus, the LHS of (B2) becomes 


|s^(Co,t)|(st(Co,i) X t) / 2'Kx^dx±_ = -27r(st(Co,t) x s^(Co,t)), 


'0 


X_L 


where we used the boundary conditions R{0) = 0,R{oo) = 1. Using this in (B2), we hnally 


obtain (10). 


Appendix C: Analytical expression of the cut-off length in terms of the vortex 
profile R{r) 


Based on the dehnition of we have: 

ln(l/^*) = fio+ lim 


f{s)ds — ln(a) 


uo 


Using equation (25) and dehnition f{s) = g{s)/s we get: 

dR{r) dR{r' 


POO /‘OO 


ln(l/^*) = /io-F lim 


'0 Jo 


dr dr' 


4rr' 


2 ^ V (r+r')'^+s‘^ 


Trio ^J{r + r'Y + s‘^ 


ds — ln(a) 


drdr', (Cl) 
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where we have used dr = 1. Notice that it is known how to compute integrals of the 

form 

/•oo K 


( 4rr' 
(r+r')^+s^ 


s ds, 


Jo A/(r + r')2 + §2 

after appropriate transformations that lead to tabulated integrals such as the ones found in 
. The integral is finite for 1/2 > z/ > 0 and has a simple pole at z/ = 0. Therefore we 


can regularise the integral over s in equation (Cl) by introducing a factor s To keep the 


integral dimensionless, we choose the following dimensionless regularising factor: (^r^) ^ • 
One can then interchange the limits a —?• oo and z/ —)■ 0, provided we can cancel the pole at 
u = 0 using a regularisation of the term ln(a). The latter can be regularised as follows: 


ln(a) = lim I In ( - I + 

^ ’ u^o+ I V 2 ' 


lya I 


—p V 


zds 


^(r + r')2 + s2 

where p is a positive number, to be chosen in order to achieve the above-mentioned cancel¬ 
lation of the pole at z/ = 0. To achieve a more familiar presentation of the integral involving 
the complete elliptic function K, we perform a standard transformation to obtain: 


ln(l/^*) = /io + lim lim 


a^oo u^0+ 


f oo POO 


'0 Jo 


dR{r) dR{r') 
dr dr' 


TT Jo 


1 +(I±^ 


K 


Arr' 2 \\ ^ 

{r + r'f^ ~ 7 (l-F)i- 


In 


r -|- r' 


—p U 


r + r' J 7(^ ^2 


ds 


dr dr'. 


We now interchange the limits, obtaining first 


ln(l/^*) = po + lim 


l/—^0+ 


fOO POO 


'0 ^0 


dR{r) dR{r') 
dr dr' 


- / K 


TT 




(r + r')2^^ ^7(1-F)i 


k-^''dk 


POO / ^ \ -pv 


zds 


'0 V' +r'J 

The integral involving the elliptic function is now in standard form, tabulated in 
Eqs.(4) and (l.a.4), with the result 


In 


dr dr'. 


y I y^ 


(C2) 


TT 


/ 4rr^ _ 2 \ 


TT 


r(7 


2 cos Z/TT r(z/ -|- 1 / 2 ) 


1 

2 F 1 V, 1; 


(r -f r'Y 


where 2 F 1 is the hypergeometric function. As for the remaining integral, it amounts to 


—p V 


+ r' J 7(^ _|_ _|_ g2 


ds = 


r(pz// 2 ) r(l /2 — pz// 2 ) 

■ 
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We see how the two integrals present simple poles at z/ = 0 coming from the Gamma 


functions. Going back to Eq.(G2), we perform a Laurent expansion on the terms in the 
square brackets, obtaining 


r^oo /*oo 


ln(l/^*) = Ho + lim / / 

^^0+ Jo Jo 


dR{r) dR{r') 
dr dr' 


P 


1 d 






In 


I 


+ 0{v) 


2’ ’’ ’ (r + r'Y 
dr dr'. 


z=0 


We see that in order to cancel the pole we need to choose p = 2. Notice that the (!2(l)-term 
that survives the limit does not depend on p, which indicates that the method is consistent. 
Therefore we obtain 


ln(l/^*) — Ho+ 


POO roo 


'0 Jo 


dR{r) dR{r') 
dr dr' 




^ y y^ 


(r + r'Y 


In 


z=0 


r + r' 


dr dr'. 


A remarkable relation allows us to simplify this: the term in square brackets above is equal 
to — ln(max(r, r')) + In2. Also, replacing po by its dehnition, equation ([ItI), we get 


ln(l/^*) — - + ln2 + 


dR{r) 

dr 


r dr 


poo /*oo 


'0 ^0 


dR{r) dR{r') ^ ^ / 

— - - - -ln(max(r, r )) dr dr . 

dr dr' 


The latter double integral can be transformed to a single integral using integration by parts. 


ln(l/^*) — - + ln2 + 


poo 

'dR{r)' 

L 

dr 


r dr + lim 

r^oo 


[R{r')f 


dr' — In r 


uo 


which is the analytical formula (12) we have aimed to prove. A numerical application 
of this formula using Mathematica gives the same result as the previous htting method, 
Eq.(26), namely = 0.3416293 ± 10“^, which validates both the numerical methods and 


the analytical formula. 
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